A Model of Correlated Fermions with dx2-y2 Superconductivity 
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Motivated by the phenomenology of the high-Tc cuprates, a two dimensional fermionic model with 
attractive interactions is here discussed. The exact solution to the two particle problem leads to 
a bound state in the d^2_y2 subspace. Numerical techniques suggest that the model has dj.2_y2 
superconductivity (SC) in the ground state at low fermionic density. Within a self-consistent RPA 
diagrammatic study, the density dependence of the critical temperature is calculated. We argue that 
in the context of d-wave SC this model fulfills the role that the attractive on-site Hubbard model has 
played for s-wave SC. We also show that another candidate, the attractive "t-U-V" model, which has 
d-wave SC at the mean-field level is actually not useful as a realization of this family of condensates 
for a variety of reasons. 
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The Hubbard model with an attractive on-site inter- 
action has played an important role in the qualitative 
understanding of s-wave superconductors. Several inter- 
esting problems, like the crossover from the BCS regime 
to the region where Copper pairs form a Bose Conden- 
sate (EC), § can be addressed within this model using 
analytical and numerical techniques. ||] Although the 
phononic degrees of freedom are not explicit, it is ex- 
pected that the qualitative properties of the attractive 
Hubbard model are the same as those of more realistic, 
and difficult to study, electron-phonon Hamiltonians. 

However, since evidence is accumulating that the high 
critical temperature (high-Tc) superconductors have a 
condensate with pairs formed in the d^^-y^ channel, 
the relevance for the cuprates of the on-site attractive 
Hubbard model is questionable. In addition, studies of 
Hc2 in cuprates suggest that the Cooper pairs are only 
a few lattice spacings in size, locating the high-Tc ma- 
terials in an intermediate region between the BCS and 
the BC limits. In this regime, the BCS mean-field (MF) 
approximation is not quantitatively accurate. |^ Thus, 
it would be desirable to have a fermionic model with a 
dx2_y2 superconducting ground state that can be studied 
with reliable non-perturbative diagrammatic and compu- 
tational techniques in the small coherence length region. 
Results obtained in this framework could be directly com- 
pared with the cuprate phenomenology. 

What model can provide the generalization of the at- 
tractive Hubbard model for the case of dx^-y2 supercon- 
ductivity? Current literature shows that many studies 
of d-wave superconductors are performed using the BCS 
MF approximation after introducing a proper attractive 
kernel in the gap equation to induce d-wave correlations. 
This approach is accurate in weak coupling where pairs 
are large, but it does not address the regime of small 
pairs which is more realistic for the cuprates. To improve 
these results, it is natural to consider as a first candidate 
for a model with d-wave SC the so-called "t-U-V" model 



where U is repulsive on-site, V is the strength of a 
density-density attraction at distance of one lattice spac- 
ing a, and t is the amplitude of a nearest-neighbors (n.n.) 
hopping term. In the MF approximation pj the phase di- 
agram at half- filling of the t-U-V model (Fig. la) indeed 
has an "island" of d-wave SC. However, note the small 
size of this phase caused by the competition with phase 
separation (PS) where electrons doubly occupy a macro- 
scopic region of the cluster to minimize the energy. A 
recent study Q has shown that the effect of PS cannot 
be simply avoided by introducing a long-range Coulomb 
repulsion since in this situation PS may be replaced by 
a charge-density-wave (CDW) state rather than by SC. 
Thus, the competition SC-PS is subtle and not created 
by the absence of long-range interactions. 

Using techniques more powerful than MF, in this pa- 
per we argue that the t-U-V model is not much useful for 
the study of d-wave superconductors. As a first indica- 
tion of this problem, let us analyze the results obtained 
using Quantum Monte Carlo (QMC) and Exact Diag- 
onalization (ED) techniques. 0] Working at small U/t 
where the MF approximation predicts a d-wave conden- 
sate, first we have studied the competition with PS. This 
phase can be observed numerically with QMC simula- 
tions since in some regions of parameter space we found 
that the mean particle density converges to two very dif- 
ferent results depending on the randomnly chosen initial 
Hubbard-Stratonovich fields. Such a behavior is typical 
of systems with two competing minima in the free en- 
ergy, as it happens in the presence of PS. In QMC the 
effect can already be seen even using small clusters of 
16 and 36 sites, and working at relatively high temper- 
atures T = t/6 to alleviate sign problems. As example, 
the presence of the two minima was numerically observed 
at U/t = 1 for \V/t\ > 0.55, reducing the MF d-wave 
stability region (Fig. la). Increasing the lattice size and 
reducing the temperature would only enhance the effect, 
since the tunneling probability would become negligible. 
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FIG. 1. (a) Mean-field phase diagram of the "t-U-V" model 
at half-filling. [5] DS,SS,PS,CDW, and SDW denote d-wave 
SC, s-wave SC, phase separation, charge-density-wave, and 
spin-density- wave order, respectively; (b) pairing cor- 

relation vs distance r obtained with QMC, on a 4 x 4 cluster, 
T = Vjl = 0.0 and V/t = -0.3 i.e. within the MF d-wave 
region. The correlations are small at distances larger than 
one. Similar negative results were obtained for other cou- 
plings inside the MF d-wave region. The error bars are smaller 
than the dots; (c) Schematic representation of model Eq.(l). 
Fermions move within the same sublattice, and interact at dis- 
tance one; (d) Binding energy Ab for the two-particle prob- 
lem solved exactly using t2o = OAtu, as example. States with 
the lowest energy in the d^2_y2, extended-s and p-subspaces 
are shown. Similar results are found for t2o < 0.5tii. If 
t2o > 0.5tii, at small |V^|/tii the actual ground state has 
p-wave, but it becomes again d^2_j,2-wave after a finite cou- 
pling is reached. 

In addition to this problem, in the regime where PS 
is not observed in our analysis, a QMC simulation still 
does not show enhanced d-wave correlations at T = t/6. 
Defining the operator that destroys a d-wave local pair as 
Af = Ci|(ci+xx -I- Ci_xi - Ci+yi - Ci-^i), with x,y unit 
vectors along the axis, then the ground state pairing 

correlation is Cd(r) = {AfAf_^_^). Fig. lb shows that 
Cd(r) decays to zero with distance. Actually, results at 
U,V = are very similar to Fig. lb. We also performed 
ED studies on 4 x 4 clusters that are in excellent agree- 
ment with the QMC data. While the QMC result does 
not rule out the stability of d-wave SC at T ^ t/6, it 
shows that (i) the small region of possible d-wave SC can- 
not be easily studied with current state-of-the-art com- 
putational studies beyond the MF approximation, and 
(ii) obviously in such a weakly correlated regime. Cooper 
pairs can only be large in size contrary to the high-Tc 
phenomenology. In addition, note also that the tentative 
d-wave SC regime of the t-U-V model occurs at half- 



filling with all electrons pairing to form the condensate. 
However, it would be desirable to have a model where a 
small density of carriers forms pairs, mimicking the ex- 
pected hole-pairing of the cuprates. The t-U-V model at 
low electronic density has s-wave SC rather than d-wave, 
Q complicating matters further. An inescapable conclu- 
sion of this analysis is that it is necessary to go beyond 
the t-U-V model for a proper study of effective fermionic 
models for d-wave SC. ||] It is remarkable that in spite 
of the recent huge effort devoted to the study of d-wave 
SC in the cuprates, the analog of the attractive Hubbard 
model for d^2_y2 pairs still seems unknown. 

The main purpose of this paper is to discuss a fermionic 
model for d^^-y'^ SC which solves the problems found in 
the t-U-V Hamiltonian. Analyzing the model studied 
here with computational techniques, it presents strong 
pairing correlations in the d-wave channel, and PS does 
not cause serious problems. The Hamiltonian contains 
an attractive n.n. density-density interaction at distance 
a, as in the t-U-V model, but it differs from it in the 
fermionic dispersion which in the new model is domi- 
nated by hopping within the same sublattice, i.e. link- 
ing next-nearest-neighbor (n.n.n.) sites. This model was 
discussed before in the context of the "Antiferromagnetic 
van Hove"(AFVH) scenario for the cuprates where the 
high Tc is induced by a large peak in the hole density 
of states (DOS) caused by antiferromagnetism. |^ The 
intra-sublattice dispersion is natural if holes move in a 
nearly antiferromagnetic background that is energetically 
costly to disturb. The attractive interaction has its ori- 
gin in AF correlations. In Ref. [|j the model was studied 
only within a MF approximation, but here we substan- 
tially improve the analysis using computational and self- 
consistent diagrammatic techniques. The Hamiltonian of 
the proposed model candidate for d-wave SC is 



H = ^eAF(k)[c[^Ck. 
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where a = A,B indicates the sublattice; eAF(k) = 
Atiicoskxcosky + 2t2o{cos2kx + cos2ky) is the dispersion; 
ill, ^20; V are parameters; and ni is the number operator, 
with the rest of the notation standard (see Fig.lc). The 
operators c satisfy anticommutation relations. They do 
not have a spin index, but carry a sublattice index which 
plays a similar role. Particles are distributed such that 
half of them are in each sublattice. Intuitively the parti- 
cles described by Eq.(l) represent "holes" in the cuprates. 

It is straightforward, and very instructive, to solve ex- 
actly the two particle problem using Eq.(l). Defining 
the binding energy as As = -E2 — 2Ei , where £"„ is the 
ground state energy of the n-particles subspace, the re- 
sults are shown in Fig. Id for extended-s, d^2_y2 and p- 
wave symmetries. The lowest energy state has A^ < 0, 
signaling the presence of a bound state, and it corre- 
sponds to c?3.2_y 2 -symmetry. At small \V\, i.e. for weakly 
bounded particles, A^ is very small in absolute value, but 
still negative. Studying the average distance between the 
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two particles we observed that at '--^ 2 in Fig. Id, 

the pair size is aheady close to its maximum value of one 
lattice spacing defining the strong coupling regime. 
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FIG. 2. (a) Two fermions in the large \V\ limit of Eq.(l) 
showing the d^2_y2-wave character of the bound state; (b) 
Two fermions in the large \V\ limit of the t-U-V model. 
The bound state is s-wave; (c) d^2_y2 pairing correlations 
vs distance r for model Eq.(l) at T=0 studied with ED tech- 
niques on a 32 site cluster. The couplings are — 1.0, 
t2o/tii = 0.4, and x is indicated. Correlations in the s— and 
p-channels are negligible. 

What is the origin of the d-wave symmetry in the 
ground state of the two-body problem? Consider the 
large \V\ limit, and let us analyze the movement of 
one particle around the other as schematically shown in 
Fig. 2a. In this regime, the energy is minimized when the 
interparticle distance is one lattice spacing at all times. 
Keeping particle B fixed at a given site, the problem now 
amounts to solving an effective four-site hopping Hamil- 
tonian of particle A moving along the four n.n. sites to B 
using the hoping amplitude tn along the diagonal. Here 
it is important to observe that the sign of tu is chosen as 
a positive number by the requirement that the minimum 
in the dispersion is at p = (7r/2,7r/2) or (0, tt), as it is 
natural in problems of holes in antiferromagnets. |^ The 
signs of and ^20 are physically relevant, unlike the sign 
of a n.n. hopping that can be changed by suitable trans- 
formations on a square lattice. Then, the ground state 
of the effective four-site problem corresponds to selecting 
a phase alternating in sign for particle A (Fig. 2a). This 
leads to a dj.-2_y-2 bound state, providing a real-space in- 
tuitive explanation for the appearance of d-wave pairs 
which complements those based on the perturbative in- 
terchange of magnons. We remark that this simple 
result found in Eq.(l) is not present in the t-U-V model. 
If the n.n.n. hopping (Fig. 2a) is replaced by the n.n. hop- 
ping of the t-U-V model, then the ground state phases of 
particle A orbiting around B at large \ V\ are as shown in 
Fig. 2b. They correspond to an s-wave bound state. 

The presence of dx2_y-2 bound states in the two body 



problem of Eq.(l) suggests SC in the same channel at 
finite particle density. However, CDW and PS states are 
also favored by a particle-particle attraction and thus 
an explicit calculation is needed to verify the existence 
of a SC condensate. For this purpose, we exactly cal- 
culated the ground-state pairing correlations C(i(r) on 
a \/32 X -\/32 site cluster using ED techniques. The 
operator is the same used before for the t-U-V model, 
simply switching the spin indices for sublattice indices. 
In order to avoid possible complications with the poten- 
tially dangerous PS regime, in the study shown in Fig. 2c 
an intrasublattice density-density repulsive interaction of 
strength \V\/^ and \V\/2 at distances ^/2a and 2a, re- 
spectively, was also included in the Hamiltonian. The 
results at = 1.0, ^20 ~ 0.4tii, and at realistic den- 

sities X = 6/32 ^ 0.19 and 8/32 — 0.25, are shown in 
Fig. 2c. The robust tail suggests strong pairing correla- 
tions in the ground state. Thus, with Hamiltonian Eq.(l) 
supplemented by mild assumptions to avoid PS it is pos- 
sible to obtain numerical signals of d-wave SC, unlike the 
results obtained before for the t-U-V model. po| , |rT| 

Hamiltonian Eq.(l) can also be studied diagrammati- 
cally using Eliashberg-type equations. |12| With this ap- 
proach, we have calculated vs x in the intermediate 
to strong coupling region. Working in Matsubara space 
and in natural units, we approximate the normal state 
proper self-energy by iteratively solving the equation 



E(k, w„) 



-E 



V;//(k-q,w„-c^:jG(q,w;), (2) 



where iV is the number of sites (we used a 32 x 32 
cluster for this calculation), T the temperature, lo^ — 
(2n + 1)ttT (—00 < n < +00), the full normal 
state one particle Green's function satisfies G(k, w„) — 
l/[zw„ - (eA_F(k) - /i) - S(k,aj„)], and V;//(k,u;„) is 
the RPA effective potential with particle-hole bubbles 
containing G rather than the noninteracting Green's 
function, to make the calculation self-consistent. Once 
the normal state G is found, for the SC state we use 



$(k,w„ 



Af(k,q, w„,a;;)$(q, uj'J, 



(3) 



where $(k,a;„) is the anomalous self-energy which 
can be considered as an order parameter for SC, 
and M{k,q,uj„,uj'J = -^V;//(k - q, w„ - uj'J 
G{q,uj'JGi-q,-uj'J. We have solved Eqs.(2,3) self- 
consistently at different temperatures and densities, us- 
ing as an energy cutoff ten times the bandwidth. The 
symmetry of the SC condensate is determined from the 
symmetry of the eigenvector of Eq.(3) corresponding to 
the largest eigenvalue, After all ring diagrams are 
summed up, M] this symmetry is d^2_y2. 
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FIG. 3. (a) Critical temperature Tc of model Eq.(l) vs 
fermionic density x, using tn = 0.165eV^, t2o ~ 0.0435el/, 
and V = — 0.075eV^, as suggested by the high-Tc phenomenol- 
ogy. [9] The solid line corresponds to the self-consistent RPA 
approximation. The dashed line is the MF or HF result; [9] 
(b) Tc vs coupling constant V, with tii,i2o fixed as in (a), 
at the density where Tc is maximum in (a). V is in units of 
Vafvh ~ — 0.075eV^, the coupling used in Ref. [9]. 

The results for Tc are shown in Figs.3a,b compared to 
the MF approximation, which in this model is equiva- 
lent to the Hartree-Fock (HF) approximation. The qual- 
itative agreement is good, and quantitatively the self- 
consistent approach reduces Tc at optimal doping by a 
factor ~ 1.5 still maintaining Tc at a high value. No 
drastic further reductions of Tc are expected by adding 
diagrams beyond RPA to the calculation. The reason 
is that in this model the vertex correction identically 
vanishes due to intrinsic features of the Hamiltonian 
Eq.(l), namely that particles move within the same sub- 
lattice, while the interaction is intcrsublatticc. To under- 
stand this effect, consider the real space representation 
of the vertex correction contribution to the self-energy 
which is given by E^'^''*'^^(r, t) c>c X^ei e2=x y ^('^ 
ei, T)G{r + ei +e2, —T)G{r + e2, t). It is clear that there 
is always a Green's function that vanishes, irrespective 
of whether r connects the same or different sublattices. 
Finally, note that the values of Tc shown in Fig. 3a are 
realistic, and the presence of an "optimal" density is a 
consequence of a peak in the DOS of e^j?(k). ^ 

Summarizing, in this paper we studied a model for 
fermions moving on a 2D square lattice with intrasub- 
lattice hopping and attractive n.n. density-density inter- 
actions. Using numerical and analytical techniques we 
conclude that in this model (i) the two-body problem 
leads to a d^2_y2-wave bound state in a natural way, and 
(ii) in the dilute limit the ground state has strong dx^-yi. 
pairing correlations. We have also provided evidence that 
the t-U-V model actually does not show a clear signal of 
d-wave SC in computational studies, and the competition 
with PS prevents the analysis of its intermediate coupling 
regime. Thus, we conclude that the new model discussed 



here is the natural generalization to d^2_y2 superconduc- 
tivity of the attractive Hubbard model. The new model 
is based on the phenomenology of the high-Tc cuprates 
which near half-filling is dominated by antiferromagnetic 
fluctuations. Phcnomenological studies of the influence 
of impurities, external fields, and other probes on dx2-y2 
superconductivity would become more accurate if model 
Eq.(l) replaces the t-U-V model. 
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